Rotational and rotationless states of weakly-bound molecules 
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Recent experimental studies of ultracold molecules 
produced by photo- and magneto-association [l]-[3] and 
of other halo molecules [4] rekindled an interest in the vi- 
brational and rotational structure of these weakly-bound 
species. In 1970, LeRoy and Bernstein [5] and, indepen- 
dently, Stwalley [6] derived a WKB "near dissociation 
equation," which related the energies of the vibrational 
levels close to threshold for dissociation to the long-range 
behavior of the potential. However, the WKB approx- 
imation is invalid for near-threshold molecular states, 
since, at threshold, the anticlassical limit is reached [7]. 
Stimulated especially by the rapid expansion of cold- 
molecule research, considerable effort has been invested 
in correcting the WKB approximation to allow for the 
treatment of weakly-bound vibrational states [8]-[15]. 

In 1972, LeRoy applied the WKB approximation to 
the rotational structure of weakly-bound molecules [16], 
which turned out to be even less accurate than for the 
vibrational levels. LeRoy concluded that "only upper 
and lower bounds could be given, rather than accurate 
predicted values for unobserved [rotational constants] 
near [the dissociation limit]." However, the rotational 
structure of weakly-bound species is of considerable in- 
terest, which is fueled mainly by the current work on 
magneto-association of ultracold atoms via higher-order 
Feshbach resonances [2], creation of ultracold molecules 
in highly-excited rotational states [17], and probing of 
halo molecules with nonresonant light [18]. 

Herein, we derive simple and accurate formulae for the 
number of rotational states supported by a weakly-bound 
vibrational level of a diatomic molecule and the rota- 
tional constants of any such levels up to the threshold, 
and provide a criterion for determining whether a given 
weakly-bound vibrational level is rotationless. 

We first introduce the quantization rule of Raab and 
Friedrich [14], [15] for rotationless states of diatomic 
molecules bound by a radial potential which behaves 
asymptotically as 
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and thus is a homogeneous function of degree —n. Our 
notation is tailored for molecules and therefore somewhat 



deviates from that of refs. [14], [15]. The quantization rule 
is given by 
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where F(Eh) is the quantization function, v is the (inte- 
gral) vibrational quantum number and Vth is the noninte- 
gral quantum number that pertains to a level exactly at 
threshold. The binding energy, = D — Ey, with D the 
dissociation energy and Ey the energy of the vibrational 
level V, is thus positive, i^b > 0. The quantization func- 
tion F{Eh) was recently derived by Raab and Friedrich 
for an arbitrary binding energy E\, [14], [15]. For a homo- 
geneous potential, it is given by: 

F(Sb) = Fth(K) + Fip(K) + Fwkb(k)] (3) 

where the individual terms, herein introduced for conve- 
nience, comprise a near-threshold dependence. 
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an "interpolation" term, which gives a smooth transition 
between 1ow-k and high-K behavior. 
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a term which corrects the reflection phase due to the 
potential of Eq. (1), 
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and a pure WKB contribution. 
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cf. Eq. (15) of ref. [14]. In Eqs. (3)-(7), the dimensionless 
wavenumber k is defined by 
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with k = \/2rnEh/h the wavenumber and m the di- 
atomic's reduced mass. The parameters b, p and u in 
eq. (3) are defined by 
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with y EE 2) (cf. Eqs. (10), (12) of rcf. [14]). Note 

that in order to avoid confusion with the rotational con- 
stant (defined below), we changed the symbol B, used in 
ref. [14], to G. The adjustable length-parameter G, which 
connects the 1ow-k and high-K behavior, was obtained in 
ref. [14] for n = 4 — 7, and is given below. 

We now turn to the case of a rotating molecule. Molec- 
ular rotation adds a repulsive centrifugal term to the at- 
tractive inverse-power potential, which gives rise to an 
effective potential 
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with J the rotational angular momentum quantum num- 
ber. For J > 0, the centrifugal term pushes the vibra- 
tional manifold due to V{r) upwards, thereby reducing 
the binding energy of the vibrational levels. For each 
vibrational level, there is a critical value, J*, of the an- 
gular momentum that pushes the level up to threshold, 
thereby causing its binding energy to vanish. Hence an 
angular momentum J in excess of J* , J > J* , dissociates 
the molecule. The effect of the centrifugal term leading 
to dissociation is shown schematically in Fig. 1. 

By making use of Eqs. (202), (204) and (205) of ref. [7], 
we can express the threshold quantum number wth for the 
effective potential, Eq. (12), as a function of J, 
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where nn(0) is the inner turning 



point, po{r) = 
lin(O) is the inner 



^y 2mVcf{{r) is the momentum, and 
reflection phase, all at threshold (zero energy). The dis- 
tance f defining the upper limit of the action integral 
must lie within a region where the WKB approximation is 
sufficiently accurate and the potential, Eq. (12), is dom- 
inated by the — C„/r" term. In such a case the sum 
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FIG. 1: (Color online) A schematic illustrating the role of the 
centrifugal term in the effective potential, Eq. (12); the en- 
ergy splittings have been exaggerated. Shown is the position 
of a rotationless vibrational level, ii(J = 0) (dashed line), as 
well as the position of the same level when pushed up by the 
centrifugal term to threshold, ^(J = J*) = Wth(J = 0) (full 
line, at threshold). When the rotational angular momentum 
J exceeds the critical value J*, the centrifugal term pushes 
the vibrational level above threshold, w(J > J*) (full line, 
above threshold), thus leading to dissociation. 



of the first two terms of Eq. (13) is independent of the 
choice of f. Since the dependence of rin(O) and (/)in(0) 
on the rotational quantum number is weak-enough to be 
neglected, the dependence of Wth on J is given solely by 
the last term. We note that 'Uth(>/) maintains the form 
given by Eq. (13) for any positive J, unlike the quan- 
tization function and other near-threshold properties of 
the potential (12), see Table 5 of ref. [7]. We now use Eq. 
(13) to derive a closed- form expression for J*, Eq. (18). 

In the limiting case, when the molecule is transferred 
from a rotationless state, J = 0, to a state with a criti- 
cal angular momentum, J = J*, the threshold quantum 
number, uth(O), decreases by: 



Az;th(J*) =«th(0)-wth(J*) 
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Since, at the same time, for a critical angular momentum 
J*, the vibrational state v is pushed up to threshold, its 
quantum number coincides with the threshold quantum 
number. 



or 



Wth(J*) = V, 
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see also Fig. 1. We thus obtain that the change of the 
threshold quantum number, Awth(</*) = Wth(O) — f, is 
nothing else than the quantization function, Eq. (3), for 
the rotationless potential V(r), Eq. (1): 



Az;th(J*) = F{E^) 



(17) 



By combining Eqs. (14) and (17) we obtain a closed-form 
expression for J*: 
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TABLE I: Comparison of the critical angular momenta J* 
and rotational constants B obtained for the three least-bound 
states of the *^Rb2 dimer in different approximations; E^, and 
B are given in 10^'' cm^^. See also Table II and text. 



V El, J* Jexact JiLB JlB B -Bexact BiLB BlB 

123 0.08 0.22 0.22 0.43 0.44 0.30 0.41 0.13 0.13 

122 97.9 4.25 4.25 4.25 4.71 4.39 4.73 4.39 3.65 

121 630.6 8.28 8.48 8.28 8.76 8.21 8.68 8.21 7.38 



TABLE II: Terms of the quantization function of Raab and 
Friedrich (RF), Eq. (3), inherent to the LB and the iLB 
approximations. Also shown are the ranges of the reduced 
wavenumber n wherein the approximations apply. See text. 
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This simple formula renders J* quite accurately, as ex- 
emplified in Table I. For instance, for the three last vi- 
brational states of ®^Rb2, Eq. (18) yields values of J* 
which are in a very good agreement with the essentially 
exact values J*xacti obtained by solving the Schrodinger 
equation with the potential energy curves taken from 
refs. [20], [21]. 

We note that Gao used the angular-momentum- 
insensitive quantum defect theory to obtain the number 
of rotational states which are supported by weakly-bound 
vibrational levels [22]. The theory, which, in general, re- 
quires evaluating the quantum defect, gives the upper 
bound on the number of rotational states supported by 
the least-bound vibrational level, Jmax = n — 2. In our 
work, Jmax is simply given by the integer part of J*, 
Jmax = Int[J*]. Hence Gao's result is seen to be a limit- 
ing case of Eq. (18) for Wth ~ v = 1. 

Neglecting any coupling of the molecular rotation, we 
can estimate the rotational constant, B, from the ro- 
tational energy, BJ*{J* + 1), required to promote the 
vibrational level bound by i?b to threshold 
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The values of the rotational constant B obtained from 
Eq. (19) for ®^Rb2 are hsted in Table I together with the 
essentially exact values, i?oxact- The latter were calcu- 
lated from 
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with the vibrational wavefunctions obtained from a nu- 
merical solution of the Schrodinger equation for the po- 
tential of refs. [20], [21]. 

From Table I one can see that for the least-bound 
state, with v = 123, whose binding energy is only 
8 X 10~^ cm~^, Eq. (19) yields a value of B which is 
by about 25% smaller than the exact one. However, for 
the two lower vibrational levels, v = 121 and v = 122, 
each of which supports several rotational states, Eq. (19) 
gives a value of B which is less than ten percent off the 
exact one. 

Table I also lists values of J* and B obtained from 
the LcRoy-Bcrnstein (LB) and the improved LcRoy- 



Bernstein (iLB) approximations. The LB approxima- 
tion results when, in the quantization function of Eq. (3), 
Fth(K) and Fcr{K) are neglected, and -Fip(K) is set to w 1. 
What remains is the WKB term which gives the semiclas- 
sical LeRoy-Bernstein (LB) quantization condition [5]. 
The iLB approximation comprises Fcr{K) and -Fwkb('*) 
but neglects Fth{K) and _F'ip(K), thereby accounting for 
short-range deviations of the true potential from V{r) of 
Eq. (1), see refs. [13], [14]. The assumptions about the 
various terms of Eq. (3) inherent to the approximations 
are listed in Table. II. 

From Table I, we see that in the case of the most 
weakly-bound state, = 123, both the LB and iLB 
approximations fail badly. The reason is the omission 
of the Fthin) term which actually determines the near- 
threshold behavior of F(£'b), see Table II. As a matter 
of fact, if for v — 123 only the -F'th('*) term were nonzero, 
we would obtain J*j^ = 0.21 and i?tii = 0.31. On the 
other hand, in the case of the v = 121 and v = 122 
levels, which are relatively far from threshold, neglecting 
the Fth{K) term is better justified. Indeed, the results of 
the iLB approximation are in good agreement with the 
accurate ones. However, since the v = 121 and v = 122 
levels are still well within the anticlassical region of bind- 
ing energies, the purely semiclassical LB approximation 
remains inaccurate. 

In the context of the present work, one can easily an- 
swer the question as to which of the weakly bound vibra- 
tional levels are rotationless, see also ref. [18]. Indeed, 
the value of the binding energy required in order for 
the first rotational state to be supported by the poten- 
tial is given by Eq. (18) with J* = 1. By combining this 
requirement with Eqs. (3)-(8), we obtain at once a cri- 
terion for a vibrational level to be rotationless, namely, 
when its binding energy satisfies the condition 



Eh < dn 
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The parameter d„ for an inverse-power potential, 
Eq. (1), depends solely on the power n. It can be ob- 
tained by solving Eq. (18) with J* = 1 numerically. The 
values of dn arc listed in Table III for ji = 4 — 7. We 
verified the validity of the criterion, Eq. (21), for the 
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TABLE III: Values of the parameter d„ appearing in the cri- 
terion, Eq. (21), which determines whether a vibrational level 
is rotationless. We also present the values from ref. [14] of the 
adjustable-length parameter G. See text. 



The values of (i„ and (i„ for different n are listed in 
Table III. 



n 


4 


5 


6 


7 


d„ 


2.8974 


1.9738 


1.6014 


1.3961 




2.8833 


1.9289 


1.5233 


1.2835 




0.73857 


0.66932 


0.63308 


0.61028 


G 


2.3528 


1.3035 


0.93323 


0.73446 



One can see that the c?J^^ parameters come close to 
the accurate ones, especially for small n, whereas the 
parameters are quite off. This is because at the binding 
energies large-enough for the molecule to support rota- 
tional states, the wavenumbcr k w 1, and thus neither 
K ^ 1 (near-threshold) nor k ^ 1 (WKB) limits apply, 
see also Table II. 



case of n = 6 by solving the Schrodinger equation for 
the Lennard-Jones (12,6) potential, and found that this 
agreed with the value of dg given in Table III within 1%. 
Also, we found that Eq. (21) provides correct predictions 
of the rotational structure of the last bound states of the 
Rb2 and KRb dimers, which we recently investigated [18] 
using the potential curves of refs. [20]-[23]. 

Within the LB and iLB approximations, the parameter 
dn can be evaluated analytically: 
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In summary, we undertook a study of the rotational 
structure of weakly-bound molecules, in which we relied 
on the quantization rule of Raab and Priedrich [14] ,[15] 
valid from the classical to the "anticlassical" limit. We 
found analytic expressions for the critical value of the an- 
gular momentum that leads to rotational predissociation; 
the number of rotational states supported by a given vi- 
brational level; the rotational constants of weakly-bound 
molecules in a given vibrational state; and a criterion for 
rotationlessness of a vibrational level. All of the above 
was found to depend just on the long-range potential 
and to check well against essentially exact calculations 
done by numerically solving Schrodinger's equation for 
the very well known potentials of Rb2 and KRb. 

We dedicate this paper to Rostislav Vcdrinskii on the 
occasion of his 70th birthday. Wc are grateful to Gerard 
Meijer for encouragement and support, and to Patrick 
Raab for discussions. 
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